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Abstract 

The redshifted 21cm line of neutral hydrogen (Hi), potentially observable at low radio frequencies (~ 50 — 
200 MHz), should be a powerful probe of the physical conditions of the inter-galactic medium during Cosmic 
Dawn and the Epoch of Reionisation (EoR). The sky-averaged Hi signal is expected to be extremely weak 
(~100 mK) in comparison to the foreground of up to 10^ K at the lowest frequencies of interest. The detection 
of such a weak signal requires an extremely stable, well characterised system and a good understanding of the 
foregrounds. Development of a nearly perfectly (-^mK accuracy) calibrated total power radiometer system 
is essential for this type of experiment. We present the BIGHORNS (Broadband Instrument for Global 
HydrOgen ReioNisation Signal) experiment which was designed and built to detect the sky-averaged Hi 
signal from the EoR at low radio frequencies. The BIGHORNS system is a mobile total power radiometer, 
which can be deployed in any remote location in order to collect radio-interference (RFI) free data. The 
system was deployed in remote, radio quiet locations in Western Australia and low RFI sky data have been 
collected.We present a description of the system, its characteristics, details of data analysis and calibration. 
We have identified multiple challenges to achieving the required measurement precision, which triggered two 
major improvements for the future system. 

Keywords: cosmology: observations - dark ages, reionization, first stars - methods: observational - meth¬ 
ods: data analysis - instrumentation: miscellaneous 


1 INTRODUCTION 

Determining what happened in the early Universe be¬ 
fore the Epoch of Reionisation (EoR) is one of the 
high priority goals of modern cosmology. A wealth of 
information on the physics of reionisation is encoded 
in the redshifted 21cm signal from neutral hydrogen 
(Hi) during these times (see e.g. IPritchard &: LoebI 
( 20121 ) for a recent review). The frequency of the 


21cm line (1420 MHz) redshifts into the low radio fre¬ 
quency range of 200 —50 MHz, for 6 < z < 30. The 
challenge of detecting the redshifted 21cm signal is be- 
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aim for a statistical detection of the Hi power spectrum 
and the Square Kilometre Array (SKA: iMellema et al 


20131) plans t o perform full tomography of the EoR 


( Wvithell2014l) . On the other hand, there are also at 


tempts to identify the EoR signature in the integrated 
spectrum of the whole sky at low radio frequencies 
- the “global” EoR signal. The physics of the cos¬ 
mological 21cm line and expectations for the global 


of authors (IShaver et al. 

I999I: Furlanetto et al. 

2006 

Morales & Wvithe 2010: 
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2010 

2012). The main expected features of the global sig- 
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nal are a trough of absorption against the CMB around 
70 MHz and a peak of emission against the CMB around 
100 MHz. Depending on the cosmological model, the 
depth of the trough is ^ 100 mK and emission peak 
~ 30 mK (Fig. [T]) against a typical sky temperatures 
of around 2000 K at ^70 MHz, decreasing to ~300K 
at ~150MHz. In the case of a sky dominated receiver 
(no noise contribution from the receiver itself), the ex- 
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Figu re 1. The expected global EoR signal jPritchard LoebI 
I 2 OI 0 I) with expected sky noise (assuming ideal noiseless re¬ 
ceiver) as a function of frequency calculated according to formula 
STsky = TskylVBr^ for frequency resolution B=lMHz, integra¬ 
tion times r = 1 h, 24h, 48 h and 96 h, and the expected sky tem¬ 
perature estimated according to formula =180K (180/i/)^ ®, 
where 1 / is frequency in MHz. 
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Figure 2. Expected signal to noise ratio derived from Figure^ 
for the same integration times. Although, the absorption trough 
at ~ 60 — 80 MHz looks like the easiest spectral feature to de¬ 
tect, the actual signal to noise ratio is not significantly lower at 
the emission “bump”. Therefore, it might be easier to detect the 
emission peak at higher frequencies, which are less affected by 
ionospheric effects and thus better understood. 


pected signal to noise ratio after 24 hours of integration 
exceeds (at the absorption trough) or is close to (at 
the emission peak) 10 at frequencies corresponding to 
the absorption trough and emission peak respectively 
(Fig. [2|). A guide of how to design an experiment to en¬ 
able data analysis capable of distinguishing the global 
21 cm signal from foreground contaminants was given 
by Liu et ahl ( 2013 1. 

In principle, measurement with a precision of the or¬ 
der of 10 mK can be achieved in ^24 hours of night 
time observations when the Galactic Centre is below 
the horizon (Fig. [T]). In practice, depending on the time 
of the year, this translates to at least a few days of 
observations to be able to average only RFI-excised 
data obt ained in optimal condit ions. S everal groups 
[EDG ES: Bowman fc Roge r^l20inl : CoRE: Chi^Eenda^ 


2009 : SAR A S: IPatra et^ l2013l S CI-HI: IVovtek et al 


2014 LEDA: iGreenhill et al.ll20l4 DARE: [Burns et al 


2012 ] are planning or operating experiments to detect 
the global EoR signal. The EDGES experiment has al¬ 
ready placed a lower limit (Az > 0.06) fo r the duration 
of reionisation ( Bowman fc Roger^l2010l l. whilst limits 
on the kinetic Sunyaev-ZePdovich (kSZ) effect from the 
South Pole Telesco pe yielded the upper limit (Az < 7.9) 
( Zahn et al ][2 oI3). 

Interest in the global EoR signal has moti¬ 
vated work to explore how well the data from 
these experiments can constrain reionisation sce¬ 
narios (Barker et al. 2012: Morandi &: Barkanall2012: 


Liu et al. 2013 : Mirocha et al.1 20131) and how instr u- 
mental effects affect the results ( Bernard! et al ][2oi3). 

A nearly perfect understanding of the signal path, re¬ 
sponse and calibration of the instrument is essential for 
the global EoR experiment and for any sophisticated 
data analysis which has to be performed to identify the 
global EoR signature. Multiple instrumental effects that 
must be controlled or under stood to high precisi o n hav e 
already been discussed by Rogers fc Bowmai] ( 2012ll . 
Perhaps the most important of these is the stability of 
the system’s response over several days of integration 
time. Integrating down to the desired precision requires 
a very stable radiometer system deployed at a very radio 
quiet site to avoid human-made radio frequency inter¬ 
ference (RFI). In order to be able to identify the global 
EoR signal in the radio spectrum the system’s response 
must be smooth and understood to a mK level, likewise 
the galactic foreground radio emission. 

Recent studies have shown that ionospheric absorp¬ 
tion and refraction can significa ntly affect the frequenc 


structure of the c osmic signal (IVedantham et al 


201 


i 


jPatta et al. 2014 1. Eurthermore, the intra-night vari¬ 
ability of the electron content in the ionosphere (even 
under “quiet” night time conditions) may impede de¬ 
tection of the global EoR spectral signature from the 
ground ( Datta et al.ll2014ll . 

This paper summarises the current status of the 
BIGHORNS experiment and discusses several top¬ 
ics (instrument, REI, data processing, calibration and 
modelling) important to every global EoR experiment. 
In Section[2]we present the BIGHORNS system, its RE- 
signal path and characteristics. In Section [3] we present 
the data collected in several remote radio quiet locations 
of Western Australia and describe processing, quality 
assessment, calibration and modelling of the these data. 
In Section |4] we present studies of the system stabil¬ 
ity and its current limitations which triggered modifi¬ 
cations to the signal path, described in Section [5j In 
Section [5] we also present the conical log spiral antenna, 
which was developed specifically for the BIGHORNS 
experiment and will be used in future deployments. 
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Figure 3. The BIGHORNS RF signal path as it was deployed in the field in 2012-2014. The antenna and the front-end (with a 3 dB 
attenuator in between) were separated by 100 m of RF-cable from the back-end electronics. The back-end comprised of a second gain 
stage (amplifiers separated by RF-filters), digitiser and spectrometer implemented on the BEDLAM board and a low power industrial 
PC for data acquisition. The three reference points marked in the image are: I - input of the front-end, L - input of the first LNA, 
O - output of the front-end and D - beginning of the digital part of the system. Unless stated otherwise, we will further refer to the 
receiver (receiver noise temperature etc.) as the entire analogue signal chain from the reference point I down to D, whilst to LNA as 
the first amplifier in the signal path (ZX60-33LN-S-I-). 


2 Instrument Description 


The BIGHORNS system is a total power radiometer. 
The main driving forces motivating the design of the 
BIGHORNS prototype were: simplicity of the signal 
path; portability; the possibility of deployment in re¬ 
mote locations without access to electricity; and low 
power consumption. The schematic of the system is pre¬ 
sented in Figure [31 The system consists of a broadband 
off-the-shelf biconical antenna (GE300E by Gompliance 
Engineering Pty Ltd) mounted 52 cm above a 3 x 3 m 
wire mesh (pitch size Ri55mm) ground screen (Fig. [4]), 
which receives radio waves from a large fraction of the 
sky. During all deployments the antenna was oriented 
in the East-West direction (Fig.|3]). 

The signal intercepted by the antenna is amplified 
in the analogue chain, which presently consists mainly 
of off-the-shelf components. The total gain of the sys¬ 
tem (~ 68 — 75 dB) was set in order to maximise oc¬ 
cupancy of the 8 bit dynamic range of the digital sys¬ 
tem with some headroom. The signal path starts with 
the front-end box. A 3 dB attenuator is inserted be¬ 
tween the antenna and the front-end box in order to 
improve the impedance match and suppress reflections 
of the LNA input noise betw een the two, as described 
by Rogers fc Bowmaiil ( 2012 ). In the case of the proto¬ 
type system with the biconical antenna, which is well 
matched only over a relatively narrow frequency band, 
the attenuator provided a better match over a wider 


frequency band. However, it is undesired because it in¬ 
creases the noise temperature of the receiver (reduces 
the signal to noise ratio) and is planned to be removed 
(improvements in the presented system will be discussed 
in Section|5l). The front-end box is powered via the Bias- 
T (Mini-Circuits part number^ ZFTB-4R2G-FT-I-) and 
houses the first low noise amplifier (ZX60-33LN-S-I- 
with gain ^22dB and noise figure (NF) ~l.ldB) and 
a two positional solid-state switch (ZX80-DR230-I-, loss 
~0.75dB). The switch is controlled by an NE555 chip- 
based circuit to switch between the antenna and a refer¬ 
ence source in order to calibrate the RF-signal in tem¬ 
perature units and also to calibrate gain variations. The 
reference source is a 5012 terminator at ambient temper¬ 
ature. The temperature is measured with a logger de¬ 
ployed in an RF-shielded box similar to the one housing 
the front-end components. The switching cadence can 
be set to two predefined configurations, Ril5.5s on an¬ 
tenna and 5s on the reference, or Ri7 s on antenna and 
5 s on reference. 

From the Bias-T the RF-signal is carried over 100 m 
of SPUMA-400-FR cable (loss ^2-7 dB in the frequency 
range 20-300 MHz) to the RF-shielded rack (Fig. |5|). A 
3dB attenuator is inserted at the end of the 100 m ca¬ 
ble in order to suppress reflected signals. Inside the rack 
the signal goes to a box where a second stage of gain 
and RF-filters are installed (Fig. |3|). The second stage 

^The codes in the brackets refer to Mini-Circuit part numbers, 
unless stated otherwise 
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of amplification consists of amplifiers ZFL-500HLN+ 
(gain ~21 dB, NF~3.9 dB) and ZKL-1R5 (gain ~41 dB, 
NF ^2.8dB). In order to band limit the signal and 
suppress high power signals below ^20 MHz, a high 
pass (41-800 MHz) filter (SHP-50+) and a low pass (0- 
300 MHz) K&L Microwave filter (Microwave part num¬ 
ber 7L120-306/T900-0/0) are used. The highest fre¬ 
quencies of interest is ^250 MHz, but having a slightly 
higher cut-off provides a potentially useful diagnostic 
band and ensures smoothness of the system’s response. 
The signal is th en passed to the BEDLAM board 
(|Brav et al.l 120131 1 , which was originally develope d for 
the LUNASKA experiment ( James et al.l l201l[ l and 
consists of analogue-to-digital converters (ADCs) and 
FPGAs configured as spectrometers. The BEDLAM 
board is a good option for a low-power-consumption 
spectrometer (< 60 W). The bandwidth of the ADC 
is 500 MHz and the signal is sampled at 960 MSam- 
ples/sec and digitised at 8 bits precision. The resulting 
sampled frequency range is more than strictly required 
given the 300 MHz lowpass filter, however the empty 
spectrum between approximately 350 and 450 MHz in 
our data provides a useful diagnostic window for vari¬ 
ations in the signal noise floor due to analogue or dig¬ 
ital effects. The spectrometer implementation uses a 
four tap polyphase filterbank (PFB) and channelises 
the data into 4096 channels with spectral resolution 
Rill7.2kHz. The power spectrum is formed and initially 
accumulated in the FPGA and the number of accumula¬ 
tions, hence time resolution, is configurable. Typically 
we have collected data with ^270 ms and 50 ms time 
resolutions. Finally, the accumulated spectra are trans¬ 
mitted over an ethernet interface to a low power indus¬ 
trial PC computer (Aleutia Ml) whi ch saves the spectr a 
on a hard drive in a FITS format ( Pence et al. 2010ll . 
The data rate at 50 ms time resolution is about 28.5 GB 
per 24 hours. The off-line data processing is described 
in Section [H 


2.1 Reflection coefficients of the components 

The scattering parameters (S-parameters) characterise 
the response of RF-components to RF-signals as a func¬ 
tion of frequency. The reflection coefficient (Sll) is a 
complex ratio of the reflected to incident voltage. There¬ 
fore, its squared magnitude (|S'llp) represents the frac¬ 
tion of power reflected and is in the range [0,1], which 
translates into (—oo, 0] in dB scale. The transmission 
coefficient (S21) is a complex ratio of voltages at the 
output and input of a device. Thus, its squared mag¬ 
nitude is a measure of gain (> OdB) or attenuation 
(< OdB). jS'llp and |5'21p are typically expressed in 
dB, and we will also typically refer to reflection coeffi¬ 
cient as magnitude expressed in dB scale. Gharacter- 
istics (such as the S-parameters or the noise figure) of 
the components were measured with a Vector Network 
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Figure 4. The BIGHORNS antenna and the front-end box (in 
the foreground) separated by 100 m of RF-cable from the camper 
trailer housing an RF-shielded rack with the back-end electronics 
(in the background) as deployed at the Wondinong Station in 
April 2014. 



Figure 5. The RF-shielded rack housing BIGHORNS back-end 
components (from the bottom to top : box with a second gain 
stage and RF-filters; BEDLAM spectrometer; power unit; and 
PC computer unit). 



Figure 6. The magnitude and phase of the reflection coefficient of 
the antenna-over-groundscreen as deployed at Wondinong Station 
measured at the reference point I (Fig. 0 . 
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Figure 7. The magnitude and phase of the reflection coefficient 
of the front-end (with a 3dB attenuator connected to its input) 
measured at the reference point I (Fig. [3jl . 



Figure 8. Measured noise temperatures of the low noise amplifier 
ZX60-33LN-S-h (blue curve) and of the entire signal path with a 
3 dB attenuator at the antenna input of the front-end with a 9^^ 
order polynomial fit superimposed (red curve). 


Analyser (VNA) whenever it was necessary. During field 
trips, the antenna’s reflection coefficient was measured 
with a portable Rohde-Schwarz ZVL VNA, whilst for 
more precise measurements (such as noise tempera¬ 
ture), which could be done in the laboratory, a Agilent 
PNA-X was used. The reflection coefficient of the an¬ 
tenna deployed at the Wondinong Station in April 2014 
(see Section IQ) was measured at the reference point I 
(Fig. [3]). Its magnitude and phase are shown in FigurelHl 
The biconical antenna is very well matched to the front- 
end (iFajp < 0.01) only in the frequency band ^172- 
184MHz and reasonably well matched (iFazp < 0.1) 
across the rest of the ~ 130 — 230 MHz band. In the 
prototype stage of the experiment the antenna was re¬ 
quired to be portable. Therefore, its characteristics were 
a compromise between what was available on the mar¬ 
ket and the predictions for the global EoR signal at the 
time of purchase. The requirements for the reflection co¬ 
efficient of the antenna will be discussed in Section 14.41 
The reflection coefficient of the front-end was mea¬ 
sured at the reference point I (Fig. [3]). Its magnitude 
and phase are shown in Figure [71 



Figure 9. The total gain of the analogue path measured with the 
Rohde-Schwarz ZVL Vector Network Analyser. A quasi-periodic 
ripple of expected periodicity 5^1.275 MHz due to a standing wave 
in the 100 m cable is shown in a frequency range 180-220 MHz. 


2.2 Receiver gain and noise temperature 

The noise temperature (Fig. [8]) and the gain of the whole 
system (Fig. |9|) were measured between the reference 
points I and D (Fig. |3|). The noise temperature of the 
system was measured in the laboratory with the PNA-X 
prior to the astronomical observations. The noise tem¬ 
perature of the terminated, standalone low noise ampli¬ 
fier is Tina ~ 80 K, the noise temperature of the entire 
signal path without a 3 dB attenuator is ~ 170 K, whilst 
the noise temperature of the entire signal path with 
a 3dB attenuator is k. 600 K. The receiver noise 
temperatures (without and with the attenuator) agree 
with the expected values calculated according to the 
Friis formula (^160 K and ~620 K respectively) for the 
system presented in Figure (3] The noise temperatures 
were parameterised by a 9*^ order polynomial in order 
to be used in the calibration process (Fig. [8]). Based 
on the reduced test, the 9*^ order polynomial was 
sufficient to describe the noise temperature data (mea¬ 
surement error ~5.4 K) in 50 — 340 MHz band. Never¬ 
theless, data analysis may be performed in the narrower 
frequency band where lower order polynomials should 
be sufficient. In a single (cold state) calibration schema, 
without the possibility of continuous determination of 
the noise temperature from hot/cold reference sources, 
higher accuracy noise temperature measurements are 
also required. 

2.3 System portability 

In order to be able to deploy the system in any remote 
location without access to electricity it was designed to 
be portable. In the field it is deployed in a modified off¬ 
road camper trailer. The whole system consumes about 
120 W of power. In order to make the system fully mo¬ 
bile, the power system consisting of a 24 V, 130 Ah bat¬ 
tery, 24 V / 23 A smart charger, and a power generator 
was assembled. The capacity of the battery is sufficient 
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to continuously run the system overnight. In the field 
the system was running on the battery during the night 
and the battery was charged for a few hours during the 
day from a power generator (Honda 20EU). 

2.4 Electromagnetic compatibility testing 

The self-generated electromagnetic interference (EMI) 
from the system was measured in an anechoic cham¬ 
ber at Curtin University. The emissions from the sys¬ 
tems were compa red against the military standard MIL- 
STD46I F (I 2 OO 7 II and commercial standard CISPR 22 
Class B (j2008ll . which define a well established proce¬ 
dure of testing electronic equipment. Moreover, the EMI 
requirements of radio-quiet sites, such as the Murchison 
Radio-astronomy Observatory (MRO) in Western Aus¬ 
tralia, are defined with respect to the military standard 
(i.e. 20dB below), and the system has to satisfy these 
requirements if it is ever intended to be deployed there. 
During the first two tests, EMI from the whole system 
was measured. In the first case the system was run¬ 
ning on the battery with the charger turned off, which 
was a typical night time operation mode in the field. 
In this configuration the system satisfies the military 
standard MIL-STD46IF ( 2007 1. However, a few nar¬ 
row band emissions in the ^ 150 — 160 MHz band and 
above 230 MHz, and wide band emissions at ^150 MHz 
and above 190 MHz were identified. Three narrow band 
emissions above 220 MHz were close the military stan¬ 
dard limit. Therefore, using these data and a model 
of radiowave pr opagat ion over the surface of the Earth 
(ITU-R P526-12 2 OI 2 II we estimated the resulting power 
that could be picked up by the antenna as deployed at 
Wondinong Station (see Section [XT]) . In this configura¬ 
tion the antenna was placed at a distance of approxi¬ 
mately 50 m from the trailer housing the RF-shielded 
rack and the power system. In the procedure we used 
sensitivity data (gain and antenna factor) measured by 
the manufacturer at 3 m distance and provided with the 
antenna. According to the antenna pattern obtained 
from the FEKO 6.3 electromagnetic simulation soft¬ 
ware, we estimated that the gain of the antenna above 
the ground screen in the direction of the rack (as de¬ 
ployed at Wondinong Station) was r; 20 — 25 dB smaller 
than the gain provided in the antenna’s datasheet. Fi¬ 
nally, we expressed the detected power in temperature 
units for 117.2 kHz frequency bins. Except for a few nar¬ 
row band spikes reaching Ri 80 — 100 mK, which should 
be excluded from data analysis, the EMI power calcu¬ 
lated at the antenna does not exceed ~10mK in the 
frequency band 50-230 MHz. Thus, we conclude that 
our measurements at Wondinong Station were not sig¬ 
nificantly affected by self-generated RFI from the rack 
during the night time observations. 

During the second test the charger was plugged into 
240 V socket and continuously charged the battery. Al¬ 


though the charger is the noisiest component of the 
system, it sat isfies the commercial standard CISPR 22 
Class B ( 2008ll . Such a configuration was not used in the 
field, where the battery was charged by the Honda gen¬ 
erator (likely the noisiest part of the system). However, 
this test allowed us to assess the EMI when the system 
is run off the 240 V socket inside a shielded building, 
which may be the case in the future. 

During the last test only the antenna, the front-end, 
and temperature probe in the metal box were inside 
the chamber and the rest of the system was outside. No 
emission could be detected in this configuration down to 
the sensitivity limit of the EMI receiver of the order of - 
10-0 dB/iV/m in the frequency band 20 — 300 MHz. It 
is very likely that these simple components are much 
quieter than the sensitivity limit. The front-end box 
is so close to the antenna that the upper limit trans¬ 
lated into temperature units in the 117.2 kHz frequency 
bin is of the order of Kelvins, which is insignificant in 
comparison to the requirements of the EoR experiment. 
Nevertheless, the measurement gave us confidence that 
there is no significant undesired emission from the “near 
antenna components”. 


3 Data collection, processing and calibration 
3.1 Data collection 

Several datasets have been collected from three differ¬ 
ent locations in Western Australia (Tab. [1}. A number 
of substantial test datasets were collected in 2012 at 
Muresk (Fig. ITUl) . While Muresk is an adequate test 
site (relatively close to Perth), the data quality is not 
good enough for global EoR studies. This is mainly due 
to RFI from relatively close transmitters in the EM 
(87-108 MHz) and digital radio and TV bands (174- 
230MHz), which can saturate the receiver or cause sig¬ 
nificant distortions in the sky power spectra. 

The quality of the data is another critical requirement 
for a measurement of the global EoR signature. There¬ 
fore, in order to collect low RFI data, several trips to 
very remote radio quiet locations were organised. Much 
better quality data, at ^270 ms time resolution, were 
collected at the Eyre Bird Observatory (EBO: Fig. ITUl) . 
The EBO is a very remote location on the coastline of 
the Southern Ocean and is hundreds of kilometres away 
from high power radio transmitters. Surprisingly how¬ 
ever, we found that even in such a remote location the 
night time data quality was affected due to long distance 
RFI propagation from distant coastline transmitters in 
South Australia, Victoria and even Tasmania (up to 
1500km away). The signals were, most likely, refracted 
via a tropospheric ducting mechanism, which allows for 
long distance propagation and is commo n in summer 


and p articularly severe near the coastline (jHitnev et al 


Il985 1. Thus, especially after dusk there were a lot of 
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Table 1 The summary of the locations and major datasets collected with the BIGHORNS system in 2012-2014 period. 


Site 

Geographic Goordinates 

Duration 

Approximate amount 
of collected data 

Muresk 

116°41’09”E , 31°44’46”S 

2012-10-04- 2012-11-05 

32 days 

2014-02-08 - 2014-02-18 

10 days 

Eyre Bird Observatory (EBO) 

126°17’52”E , 32°15’08”S 

2013-12-07- 2013-12-17 

7 days 

Wondinong Station 

118°26’24”E , 27°51T0”S 

2014-04-04- 2014-04-11 

5 days 



Figure 10. Outline map of Western Australia with Perth and 
the three deployment sites (Tab. ^ marked. 


RFI signals in the FM and digital TV bandsH Par¬ 
ticularly, the later were undesired as they affected the 
most promising band, where the biconical antenna is 
best matched. Hence, the third site was inland at the 
Wondinong Station (Fig. [TOl), where the best data to 
date were collected. 

3.2 RFI excision and data reduction 

The main requirement for the BIGHORNS deployment 
site was for a low RFI environment. Thus, the col¬ 
lected data are intrinsically relatively clean, even within 
the FM band. The main source of the high power 

^We found some tropospheric ducting forecast services to be 
good predictors of the level of RFI seen in the data. E.g. 
http: //www. dxinf ocentre. coin/tropo_aus .html 


Total power in 0-480 MHz band 
Threshold value (= median + 0.0002 mW) 



Local time 

Figure 11. Total power in the 0 — 480 MHz frequency band ob¬ 
served over 12 hours at the Wondinong Station with a cut-off 
threshold value. 


RFI signals are from the low orbit telecommunication 
ORBCOMM satellites at around 137.5 MHz (right-hand 
circularly polarised) and from airplanes in the frequency 
band 117.975 — 136 MHz. The power transmitted by 
these devices is relatively high and can saturate the 
receiver. Very high power in a narrow band can also 
cause distortions in the noise floor of the remaining 
spectrum. In order to efficiently exclude RFI affected 
data they should be collected at the highest possible 
time resolution, which is limited by the maximum data 
transmission and storage rates in the system. Therefore, 
the Wondinong Station data were collected with 50 ms 
resolution. 

Due to the fact that the calibration switch oper¬ 
ates autonomously, the first data processing step identi¬ 
fies integrations collected on the antenna and reference 
source and saves this information to the FITS file head¬ 
ers. The state identification procedure calculates total 
power in every integration in a certain frequency band 
(60— 110MHz for the Wondinong data), compares it 
with a threshold value and flags as an antenna or refer¬ 
ence integration. The threshold value is typically calcu¬ 
lated automatically (based on the distribution of total 
power) or can be pre-defined according to observed to¬ 
tal power for each of the states. When a state transition 
is detected, three surrounding integrations are flagged 
as undefined in order to be excluded from further anal¬ 
ysis. In the next processing step, in order to excise data 
affected by RFI, the following criteria are applied: 
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percentage of excised data 
power in a singie channei > -35 dBm (criterion 1) 
broadband excisiori (criterion 2) 


AOfiagger software 


MB 1: - : - : - : -I 


(criterion 3) 





100 
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250 300 
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Figure 12. Percentage of all the Wondinong Station data excised by the three criteria in the frequency resolution of 117.2 kHz. The 
excision rates of each individual criteria were calculated independently of the other criteria. The power based criteria are constant 
across the frequency at 5^13% (red curve) and 5^16% (blue curve). The heavily RFI affected channels were almost entirely rejected by 
the AO-fiagger software. Therefore, the green curve lies on top of the black curve in these channels. 





150 200 250 300 


Frequency [MHz] 


Figure 13. Comparison of RFI occupancy between the Wondinong Station, the EBO and Muresk. Please note that RFI threshold 
values used in criteria 1 and 2 (see text) for the EBO and Wondinong data would reject all the Muresk data. Thus, the figure shows 
occupancy obtained after cuts appropriate for the Muresk data were applied, which situates the “rejection floor” slightly below the 
other two, but also shows a lot of narrow band emission. 


1. Single channel excision - if the power measured 
by the spectrometer exceeded -35 dBm in any sin¬ 
gle channel the whole integration was excised. The 
value of the threshold was determined empirically 
in the laboratory with a signal generator connected 


to the antenna input of the front-end. It has been 
observed that power >-35 dBm affects the noise 
floor by >1%. The threshold was set at -15 dBm 
for the Muresk test data. The average excision rate 
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0.00111 
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0.00012 


Figure 14. Prediction for sky model brightness temperature Tu{0, (p) with over-plotted contours of simulated antenna pattern Pu(0, (p) 
(left panel) and their product (right panel) at z/=150MHz. The upper pair of images corresponds to Galactic Centre below the horizon 
(at 2014-04-06 10:34:51 UT, corresponding to approximate local sidereal time 07^27^), and the lower images correspond to Galactic 
Centre transit (at 2014-04-06 21:07:51 UT, corresonding to approximate local sidereal time 18^02^). Summation of values in all the 
pixels of the original images (before the orthographic projection) in the right panel leads to a MHz) as in the corresponding 

curve in Figure fT^ 


of this criterioid for the data from Wondinong Sta¬ 
tion was ~13%. 

2. Broadband excision - if the sum of power in 
the entire band was above the median galactic 
noise level plus the threshold of -37.5 dBm (corre¬ 
sponding to 5 X 10® in the original arbitrary power 
units in which the cut is performed at the early 
data processing stage) the integration was excised 
iFig. fTTll . For the Muresk data the threshold value 
was -14.5dBm (10^^ in arbitrary units). The pur¬ 
pose of this criterion is to excise integrations af¬ 
fected by high RFI power appearing in more than 
one channel. The average excision rate of this cri- 


^The average excision rates of every criteria were calculated in¬ 
dependently of the other criteria 
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terion for the data from Wondinong Station was 
^16% (Fig. 112)). 

3. AO-flagger - all t he data were flagged with the 
AOflagger software (lOffringa et al.ll2012l) . and the 
flagged data were excluded from further process¬ 
ing. The “rejection floor” of this criterion is ~ 
1.3 — 2% in the frequency band 60 — 300 MHz for 
the data from Wondinong Station and reaches 
100% for RFI affected channels iFig. [T2|) . 


The average excision rates of each criterion applied 
to the data collected at Wondinong Station are shown 
in Figure fT2l The RFI occupancy comparing the three 
sites, where the system was deployed in years 2012-2014, 
is shown in Figure [T21 After the RFI excision process, 
the data volume was reduced to a more manageable vol- 
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ume by averaging every N integrations on the antenna 
and a reference load, where N was 700 for the sky signal 
and varied between 200 and 300 on the calibration sig¬ 
nal. The number of averaged single integrations in every 
frequency channel was saved to a separate FITS file for 
later analysis. Further data analysis was performed on 
the reduced dataset. 

3.3 Sky signal modelling 



In orde r to calculate the expected sky spectrum, a sky 
model (|de Oliveira-Costa et al.l 1200811 was integrated 
with a simulated antenna pattern. A model of the bi- 
conical antenna 52 cm above a 3 x 3 m ground screen 
and soil was developed in the FEKO 6.3 package in or¬ 
der to calculate the expected antenna pattern, Pu{0, (l>), 
and efficiency, 'q{v). The coordinates are angles 

in the antenna frame and can be directly translated to 
horizontal coordinates (az,el) according to orientation 
of the antenna in the field. For a given time and fre¬ 
quency, v, the antenna pattern P^(9,4i) was integrated 
over the entire sky with the sky brightness temperature 
T^{9, (/)), obtained from the sky model, according to the 
formula: 


Figure 15. Efficiency ( 77 ) of the biconical antenna obtained from 
the FEKO model. The dip at P:il30MHz is caused by ground 
screen related effects, and it is also observed in the data, although 
slightly shifted in frequency (Fig. fl^ . Using the FEKO model it 
was verified that the depth of dip depends on the height of the 
antenna above the ground screen and that the dip is not present 
when the antenna is in free space (without the ground screen). 


spectra observed when the switch is set to antenna and 
load can be expressed by the following equations: 


Pref=9{Tamb{l-\Tl\^)+T;,,) (3) 


!i^PA9A)d^ 


( 1 ) 


The above formula is time dependent through the 
implicit time dependence of T^{9,(j)). Many equations 
and quantities in this paper are frequency and time de¬ 
pendent, but their frequency and time dependence is 
usually not explicitly shown for brevity. 

Predictions for the sky model, antenna pattern, and 
their product in zenith orthographic projection are 
shown in Figure [TH 

Resistive losses in the antenna rods and due to signal 
picked up from the ground were estimated in the FEKO 
simulation (Fig. [121). Thus, the expected output signal 
from the antenna (before the balun) is 


Pant — QPrec ( 4 ) 

where Pref and Pant are the powers measured on the 
antenna and reference respectively, and g is the system 
gain post low noise amplifier (LNA) at the receiver. 

T^av is receiver noise temperature when the RF- 
switch is set to the 5017 terminator measured with the 
VNA between points L and D (Fig. [3]), and [F/p is the 
reflection coefficient between the receiver and reference 
source (in this case a 5017 resistor). 

Tree is the sum of signal collected by the an¬ 
tenna {Tant) and transmitted to the receiver (not re¬ 
flected back) and full (including noise reflections in the 
antenna-receiver system) receiver noise (Trevn)' 


Tjnodeli^) — Tpatterniy^'n^^^ “b (4 V{^')')Pamb: (2) 

where Tamb is ambient temperature. Further losses in 
the signal chain (before the RF-switch inside the front- 
end), due to balun and cable inside the antenna, were 
included in the calibration process. A comparison of cal¬ 
ibrated and modelled spectra of the sky at Wondinong 
Station when the Galactic Centre was transiting and 
when it was below the horizon are shown in Figure [T6l 


3.4 Absolute calibration of the signal 

In the existing calibration schema a two position switch 
was used, switching between the antenna and a 5017 
resistor at ambient temperature (cold load). The power 
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Tree — Tant{^ — \^al\‘^)L + Trevn, (5) 

where L is a frequency dependent loss between the an¬ 
tenna terminals and the first LNA (due to losses in the 
RF-switch and an optional attenuator). Losses of these 
components as a function of frequency were measured 
with the VNA. IFq/P is the reflection coefficient between 
the antenna and receiver. Hence, the term (1 — iFajp) is 
the fraction of antenna signal entering the receiver (not 
reflected back to the antenna). Tai can be expressed in 
terms of impedances of the antenna {Za) and receiver 
(Zi) as: 


Za + Zi 


( 6 ) 
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Tant in equation[5]is the antenna temperature (in loss¬ 
less antenna case Tant ~ Tsky)- Because frequency de- 
pendance of losses inside the off-the-shelf biconical an¬ 
tenna could only be estimated according to datasheets 
of the balun and cable, non-zero losses inside the an¬ 
tenna {La ^ 0) were optional in the calibration process 
as they might be a source of significant uncertainty. Af¬ 
ter taking into account non-zero losses inside the an¬ 
tenna, Tant can be expressed as: 


Tant — TskyLa -f (1 — La)Tamb- (7) 

Finally, the receiver noise {Trcvn) in equation [5] can 
be expressed as a sum of the standard El receiver noise 
temperature (measured between points I and D 
in Figure [3|), and a reflected noise term {Tnw) due to 
receiver input noise reflected back and forth between 
the antenna and receiver (“noise waves”) resulting from 
the impedance mismatch of these elements (Fig. [3]): 


T i _ rpa 

rc-vr). — -L 


Tn 


( 8 ) 


In the first approximation Tj-cvn ~ T^^v^ bnt higher 
accuracy calibration requires taking into account re¬ 
flected noise terms {Tnw)- The biconical antenna 52 cm 
above the ground screen (Section[2]) is very well matched 
to the front-end (iFaip < 0.01) only in the frequency 
band ^172-184 MHz. Thus, the contribution of reflected 
noise is substantial outside this range. The contribu¬ 
tion of the reflected noise was estimated with an ap¬ 
proach described bv IHogers fc Bowman (2012), where 
it is treated as noise waves due to LNA input noise re¬ 
flected back and forth between the antenna and LNA. 
The noise wave amplitudes were suppressed by a 3 dB 
attenuator inserted between the antenna and input of 
the front-end. However, this is undesirable because it 
significantly increases the noise temperature of the re¬ 
ceiver (improvements in the presented system will be 
discussed in Section [5|). The procedure of determining 
noise wave contribution was described in detail in the 
aforementioned paper. In summary, the amplitudes of 
the noise wave terms were estimated by measuring a 
power spectrum of an open (or shorted) cable of suffi¬ 
cient length (in our case ^ 10 m) connected to the input 
of the front-end (instead of the antenna); and measuring 
the reflection coefficient of this cable; and the reflection 
coefficient of the front-end (which was measured any¬ 
way) . The observed power spectrum of the cable can be 
described by equation [51 where Tant is set to Tamb (as 
the cable “sees” the physical temperature), and Trcvn 
can be expressed in terms of correlated (Tc, Tg) and un¬ 
correlated (T„) portions of the noise from the input of 


^We refer to standard receiver noise temperature as measured 
with the VNA which uses a well matched reference source. In 
contrast to real, mismatched sources such as an antenna. 


the LNA as in equation 8 in Rogers fc Bowm^ ( 20I2h . 
After assuming certain parameterisation (e.g. linear) of 
the Tc,Ts and T^ terms, their parameters can be fitted 
in order to reproduce the observed power spectrum of 
the cable. 

An alternative, anechoic chamber based method of 
determining full receiver noise (including reflected noise 
contribution) was also developed and tested. It is based 
on the assumption that the brightness temperature in 
the shielded (~80 dB at the frequencies of interest) ane¬ 
choic chamber is equal to the ambient temperature. A 
proper analysis of how well the anechoic chamber ap¬ 
proximates a black-body is considered as very impor¬ 
tant if the method is going to be further used and im¬ 
proved, but is beyond the scope of this paper. How¬ 
ever, a first order approximation can be obtained based 
on normally incident plane wave and reflectivity (de¬ 
viation from a perfect absorber) of the absorber (FS- 
600H), which is ~I0% at 20MHz, ~3% at 30MHz and 
^1% above 80 MHz. Thus, we estimate that the cham¬ 
ber acts as a black-body within a few percent down to 
^30 MHz and even better (within ^1%) at frequencies 
above 80 MHz. Equation [5] can then be converted to: 


= Tree — Tamb{^ — \^al\‘^)L, 


(9) 


where all the values on the right hand side can be 
determined: Tree according to equation [151 Tamb with 
a temperature probe; and coefficient Fq; can be cal¬ 
culated from equation [51 after measuring the reflection 
coefficient of the antenna in the anechoic chamber. In 
order to measure Trcvn the system was deployed in the 
anechoic chamber in almost exactly the same form as in 
the field. The resulting values of Trcvn were later used 
to calibrate the sky data. T he method yields very sim - 
ilar results to the method of [Rogers &: Bowman! (l2012i) 
and might be further improved and used in the future. 

Equations [4][8| can be solved to determine Tree which 
can be expressed as 


T — 

rec — 


Pa- 


P 


ref 


-{Tamb{l-\ri\^)+T^cv)- ( 10 ) 


The factor (1 — |F;p) can be approximated by I, be¬ 
cause a 50II terminator is very well matched to the 
front-end input. Finally, the calibrated sky tempera¬ 
ture, corrected for reflections in the antenna-LNA sys¬ 
tem, can be obtained from equations [SJ [3 [5] and [TUI as 


rr _ Trcvn ^ 1 

{l-\Val\^)LLa^^ lJ ^ > 
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3.5 Statistical errors 

The statistical error on Tree can be expressed by the 
following formula: 



( 12 ) 


where STamb is the error on the temperature mea¬ 
surement and SPref and SPant are statistical errors 
of the power measured on the reference and antenna, 
respectively. Equation [T^] holds under the assump¬ 
tion that the errors {SPref, SPant and STamb) are in¬ 
dependent and Gaussian distributed, which we veri- 
hed on subsets of the data. It can be simplified un¬ 
der the assumption that the error on the measured 
power satisfies the radiometer equation SP/P = Bt, 
where B is the resolution of the frequency bin and r 
is the integration time. If a long integration time r is 
achieved by averaging M single integrations of duration 
Ti then the error on the temperature measurement also 
decreases as STamb = cramblVM = Uamb /where 
Oamb ~0.03K is the estimated error on a single mea¬ 
surement with the presently used temperature probe. 
Consequently, the error on such an averaged Tree can 
be expressed by 


STree = 


Pa 




amb 


P 


ref 


B 


^'^ant ^Tef 


' amb 


(13) 

where Tant and Tref are integration times on antenna 
and reference, respectively (in a situation when they are 
not the same). 


3.6 Comparison of calibrated data and sky 
model predictions 

An example of two calibrated integrations (average of 
700 single 50 ms integrations) when the Galaxy was 
transiting and below the horizon is shown in Figure ITHl 
The relative difference between calibrated spectrum 
and the sky model is below 10% in the 70 — 200 MHz 
band, and gets down to 5% when an extra 0.4 dB 
of unaccounted attenuation is added to the attenua¬ 
tion inside the antenna. These values are within the 
sky model error, which was e stimated by its authors 
( de Oliveira-Gosta et akllioO^i to be Ril0% at low fre¬ 
quencies (20 — 400 MHz). 

The discrepancy may be attributed to an overesti¬ 
mation of the antenna efficiency ( 77 ), underestimation 
of the signal loss between the antenna and the LNA, 
incorrect subtraction of reflected noise contribution, or 
uncertainty in the sky model itself. In order to test our 

PASA (2015) 

doi: 10.1017/pas. 2015.XXX 



Figure 16. Calibrated data collected on 2014-04-06 at the 
Wondinong Station (with RFI affected channels removed) com¬ 
pared with the sky model integrated with antenna pattern. The 
dash-dotted lines show a calibration with an extra 0.4 dB loss in¬ 
side the antenna, which could be the case due to lack of exact 
characteristics of the antenna. 


understanding of the signal calibration without relying 
on the exact knowledge of the additive noise terms Trevn 
(mainly uncertainty of reflected noise), the difference of 
the two calibrated integrations was compared with a 
difference of the corresponding sky model predictions 
(Fig. ITTI) . Figures ITHl and ITTI show that an accurate ab¬ 
solute calibration requires a very good knowledge of the 
antenna efficiency, which is a poorly constrained quan¬ 
tity. In the presented analysis the value of 77 was ob¬ 
tained from the FEKO simulation, and in practise it is 
difficult to have a better estimation. The antenna ef¬ 
ficiency is typically measured with a ^ 20 % precision, 
5-10% is considered good and the best accuracy tha t 
can be achieved is ~1% ( Johnston fc McRorv 1998h . 
Even such an extremally good measurement would in¬ 
troduce a frequency dependent error >1 K outside ^ 
130 — 170 MHz band (at ~150MHz “cold sky” temper¬ 
ature equals typical ambient temperature thus antenna 
efficiency does not matter). Although it might not be 
possible to absolutely calibrate the signal to a mK pre¬ 
cision, as long as the antenna efficiency is a smooth 
and “well behaved” function (without bumps and wig¬ 
gles) of frequency, the global EoR signature can still be 
detected. The precise absolute calibration also necessi¬ 
tates good characteristics of losses inside the antenna 
which, in the case of the off-the-shelf biconical antenna, 
could not be measured precisely (without breaking the 
antenna) and only approximate data-sheet values were 
used. 


3.7 Data quality assessment 

The calibration schema was used to calibrate 5 days of 
averaged antenna integrations by the corresponding av¬ 
eraged reference integrations and saving the resulting 
dynamic spectra to FITS files. In order to assess the 
quality of the antenna data, a median dynamic spec- 
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Figure 17. The difference of the calibrated spectrum at two dif¬ 
ferent epochs (21:07:51 and 10:34:51 2014-04-06 UT) compared 
with the corresponding differences of the sky model. The red dash- 
dotted line corresponds to difference of calibrated data with an 
extra 0.4 dB (Fig.^. The two model curves were generated for 
antenna efficiency rj=l and r] as from the FEKO simulation. 


trum for 1 sidereal day was generated by grouping all 
data into local sidereal time (LST) bins of size approxi¬ 
mately 1 min and taking the median of each bin. Every 
integration of the entire 5 day dataset was then divided 
by the corresponding value (with the closest LST) of the 
median dynamic spectrum. An example piece of the nor¬ 
malised dynamic spectrum is shown in Figure 1181 The 
daytime data are significantly affected by solar activ¬ 
ity, but the night time data are stable (typically to less 
than ^1%). Several RFI related features marked in the 
image are due to: military satellites at ^235-290 MHz, 
FM band at 87.5-108 MHz and ORBCOMM satellites 
at 137-138 MHz. The relatively high power transmitted 
by ORBCOMM satellites sometimes caused some unde¬ 
sired leakage-like effects on the nearby channels (marked 
with arrows in the bottom of the plot), which were not 
excised by the earlier RFI flagging criteria, but were 
identified and excluded based on the normalised dy¬ 
namic spectrum. The first 24 hours were also affected by 
thunderstorms within some 100 km. The ft!l.275MHz 
ripple resulting from the reflections in the 100 m cable 
is present almost in the entire band, but it is well visi¬ 
ble as vertical lines below 100 MHz due to a particularly 
bad match at these frequencies. Some horizontal lines 
are due to ORBCOMM power affecting the entire band 
which may enforce an even more rigorous threshold in 
the RFI excision criteria 2. The normalised dynamic 
spectrum is a very useful diagnostic tool which enables 
identification and excision from further analysis integra¬ 
tions still affected by RFI or other undesired effects. 

4 Limitations of the system 

Averaging of the spectrum to mK precision requires col¬ 
lecting data over several days, as only a few hours of 
the best data per day can be used. Such a procedure 
requires a stable system. Stability of the BIGHORNS 
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Figure 18. Normalised dynamic spectrum of approximately 
66 hours. Frequency scale is on horizontal axis and time flows 
from the bottom to the top. The first night was affected by thun¬ 
derstorms within some 100 km. Nevertheless, excluding frequency 
bands affected by RFI (marked in the figure) due to FM (87.5- 
108 MHz) and ORBCOMM (137-138 MHz and sometimes nearby 
channels), the night time data is mostly stable to <1% in the fre¬ 
quency band 50-235 MHz. Several other features such as: daytime 
solar activity and RFI due to military satellites (~235-290 MHz) 
were also marked in the image. The P=il.275MHz ripple is well 
visible as vertical lines at frequencies below 100 MHz. 


system was tested in the laboratory and in the anechoic 
chamber. 

4.1 Stability of the reference data 

During the laboratory tests the antenna input of the 
front-end was terminated with a 50D resistor. A large 
laboratory dataset (475 hours) was reduced by averag¬ 
ing all 200 integrations (^270 ms each) in every single 
FITS file (^54.6 sec) and saving them to a separate ont- 
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Figure 19. Fit of constant temperature in the 155-190 MHz band 
to calibrated reference integration of p^ 54.7 hours duration (lower 
panel), and a distribution of residuals with a Gaussian fit su¬ 
perimposed (upper panel). The standard deviation of the fitted 
Gaussian represents ST^efi'^ ~54.7 hours) ^ 3.8 mK in Figure f2Q] 


Figure 20. The error (r) in a single frequency bin of 

117.2 kHz within 155-190 MHz frequency band (derived as de¬ 
scribed in the text) as a function of integration time t. Every 
point in this figure was obtained from the Gaussian fit to a dis¬ 
tribution of residuals IFig. I19|l . 


put FITS file. The variations of the recorded signal were 
estimated by normalising the averaged dynamic spec¬ 
trum to the hrst or a median integration. The response 
of the system changes by a few percent due to tem¬ 
perature variations and changes in the system’s gain as 
a function of temperature. Switching between the an¬ 
tenna and a reference source allows us to calculate the 
instantaneous gain and thus eliminate the effects of gain 
variations from the calibrated data. 

Using the same large laboratory dataset, the ability 
of the system to integrate down according to the ra¬ 
diometer equation was tested. Because spectra collected 
on the terminated antenna input are indistinguishable 
from those collected on the reference, odd integrations 
were calibrated with the even ones. Such calibrated data 
were later averaged and a constant room temperature 
was fitted in the 155-190 MHz sub-band to the resulting 
average spectrum. The Gaussian distribution of residu¬ 
als of this fit yielded a statistical error STref ~1.7mK 
of the mean integration resulting from averaging all 
the integrations, which is within 0.1 inK of what is ex¬ 
pected from equation [12] and r Ri (475/2) hours and 
B=117.2kHz. The same procedure of fitting Gaussian 
to a distribution of residuals (Fig. [191) to determina¬ 
tion of the error as a function of integration time which 
is shown in Figure |20l 

The stability of the data collected at Wondinong Sta¬ 
tion was also assessed in the same way. Every 200 refer¬ 
ence integrations were averaged to a single (~10sec) 
integration, and odd integrations were calibrated by 
the even ones. All the calibrated reference integrations 
(Ril3.2hours) were averaged and a constant (across fre¬ 
quency) temperature ht was subtracted. The error of 
Ri 7.4 mK was derived from the distribution of resid¬ 
uals and agrees very well with a value ~ 7.3 mK ob¬ 
tained from equation [13] for a given integration time 
and {Tamb + T//^)=386K. The dependence on the inte¬ 


gration time is also in a good agreement with the ex¬ 
pectations. 

Finally, the quality of the reference data was always 
verified by normalising the averaged dynamic spectrum 
of the reduced reference data to the first (or median) 
integration, which would allow identification of any un¬ 
desired variations in the reference signal and system’s 
response. 

4.2 Stability of the antenna data 

In order to test the stability of the antenna data, the 
system was deployed in the anechoic chamber in almost 
the same form as in the held, and a few days worth 
of data were collected. The uncalibrated antenna and 
reference power varies by ^1%. The calibrated signal 
varied by ^0.1% which corresponds to variations of the 
physical temperature in the chamber. The data were 
collected at ~270 ms resolution and were reduced as de¬ 
scribed in Section [3l Single integrations were averaged 
resulting in a reduced dynamic spectrum at ^54.6 sec 
resolution (200 single integrations) and reference at 
~15.8sec (43-65 single integrations) resolution. The re¬ 
sulting dynamic spectrum was calibrated to Tree accord¬ 
ing to equation [5] Then, for every frequency channel a 
standard deviation of Tree was calculated from a sample 
of increasing number of integrations (Fig. I21|) . 

The standard deviation converges to a constant value 
after a sufficient number of integrations and remains 
constant at least over a period of ^5 days of the test 
duration. Its value, representing observed STrec of a sin¬ 
gle frequency channel of a single calibrated mean inte¬ 
gration {Ti R!54.6sec), varies across frequency between 
0.25-0.4 K. The expected values of STrec were calcu¬ 
lated according to equation [131 and a good agreement 
between the observed and expected values of STrec is 
shown in Fig. [22] The good understanding of the statis¬ 
tical error on the calibrated data collected in the ane- 
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Figure 21. Standard deviation of the Tree as a function of num¬ 
ber of sample integrations used to calculate it for selected fre¬ 
quency channels. 



Figure 22. Standard deviation of Tree calculated for all integra¬ 
tions of the chamber data fFig. f^Tll compared with the expected 
value of STree calculated according to equation 1131 


choic chamber enabled us to estimate expected statisti¬ 
cal errors for the sky data. 

The expected STrec for the sky data was calculated 
according to equation |T3| where the ratio Pant/Pref for 
a “cold sky” was taken from the Wondinong data. The 
values of STrec for single 117.2 kHz channels calculated 
for several integration times (including a single aver¬ 
age integration of ~35sec) are shown in Fig. |23| The 
expected error of 10 average channels, corresponding 
to a frequency resolution of ~ 1.17MHz, would be \/T0 
smaller. 


4.3 Limitations of calibration 

The two positional calibration schema works under the 
assumption that noise and S-parameters characteristics 
of the LNA do not vary significantly with variations of 
environmental parameters. Modern LNAs satisfy this 
assumption fairly well. The LNA response was tested 
in the laboratory by heating up (to ~320-330K) and 
cooling down (to ^276-283 K) the first amplifier in the 
chain and observing changes in measured noise tem¬ 
perature and S-parameters. The maximum differences 
of the LNA’s parameters between cold and hot states 
were: gain RiO.OSdB; reflection coefficient RiO.SdB; and 


Figure 23. The expected statistical error STree calculated ac¬ 
cording to equation 1131 as described in the text for the case of 
“cold sky” (when the Galactic Centre was in the lowest point 
below the horizon). The error does not significantly increase at 
lower frequencies, as could be expected from the high sky tem¬ 
peratures (Fig. due to the very poor antenna match at these 
frequencies (refiecting back more than half of the sky signal). 


noise temperature ~1-2K. The parameter changes were 
affected by a systematic shift (without any significant 
frequency structure), which could result in a systematic 
error of the absolute calibration ft!l-2K. Daily varia¬ 
tions of temperature at Wondinong Station were ^288- 
316 K over the period of data collection, and only ^294- 
306 K over 4 hours after the sunset (before the rise of 
the Galactic Centre). Thus, the error on the absolute 
calibration was likely smaller than the estimations from 
the laboratory tests. Particularly, for the best parts of 
the data (a few hours after sunset) when the temper¬ 
ature was close to room temperature at which LNA’s 
characteristics was measured in the laboratory. A new 
front-end equipped with an in-built “hot and cold” ref¬ 
erence source will enable us to continuously determine 
the instantaneous receiver noise lSection l5.ll) . 

The main difficulty in calibrating the data collected 
with the presented system comes from the effects of 
source impedance (antenna or 50 H terminator) propa¬ 
gating further downstream the signal path. Particularly, 
the aforementioned Ril.275MHz quasi-periodic ripple 
due to reflections of the signal on both ends of the 100 m 
cable connecting the front-end with second gain stage 
cannot be fully calibrated out and can only be mod¬ 
elled. This is due to the fact that it has a different am¬ 
plitude and phase as a function of frequency when the 
RF-switch is set to the antenna and reference. These 
effects can be diagnosed and quantified by measuring 
the reflection coefficient at the output of the front-end 
(S22) at the reference point O (Fig. [3]) which is signif¬ 
icantly different when the RF-switch is set to antenna 
and a 5011 terminator (Fig. [24l) . Several modifications 
in the front-end, which should efficiently eliminate this 
problem, have already been tested and will be discussed 
in Section EH 
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Figure 24. The reflection coefficient (S22) of the front-end mea¬ 
sured at point O (Fig. |3j with the three different sources con¬ 
nected to its input. 


4.4 Antenna mismatch 

The off-the-shelf biconical antenna, as presented in 
Section 12.11 is very well matched (reflection coeffi¬ 
cient ^ -20 dB) over a relatively narrow frequency band 
(~172-184MHz). Nonetheless, very good match in a 
wide frequency band (~50-100MHz) facilitates the cal¬ 
ibration procedure and thus is a very important fac¬ 
tor for the global EoR experiment. A better matched 
antenna together with the “isolation” in the front-end 
(Section l4.3ll can further reduce the downstream propa¬ 
gation effects. While the effects of the LNA input noise 
reflected between the antenna and LNA cannot be en¬ 
tirely eliminated, they can be significantly reduced with 
a very well matched antenna. Assuming that all the mis¬ 
match effects downstream from the hrst LNA can be 
calibrated out by cold or hot-cold reference (discussed 
in Section 14.3p , some requirements for optimal reflec¬ 
tion coefficients come from the considerations of noise 
reflected in the antenna-receiver system and absolute 
calibration (Section 13.4|) . In the first approximation 
(single reflection of the receiver input noise from the 
antenna) the amplitude of these reflections is propor¬ 
tional to the product of the antenna reflection coefficient 
(ILop) and input noise temperature of the receiver (typ¬ 
ical value of the order of Tin ~50K). The requirement 
for ri„|rap<lmK leads to ITqP <-47dBm - nearly a 
perfect match. Such a good match cannot be achieved in 
practise over a finite bandwidth. In the antenna engi¬ 
neering a <-10 dB match is considered good, <-15 dB 
match very good and <-20 dB in a wide frequency 
band (e.g. 100 MHz) is extremally good and difficult 
to achieve. Usually a <-30dB match can be achieved in 
a narrow frequency band of the order of 5-10 MHz (e.g. 
biconical antenna in Figure [H). The condition might be 
relaxed by requiring only smoothness of the reflection 
coefficient in a restricted frequency range (tens of MHz) 
to satisfy the condition Tin\^Ta\'^<l mK, which results 
in condition IATqP < 2 x 10“®. This means that the 
smoothness of the reflection coefficient itself should be 


smaller than 10“^, which is also difficult to achieve 
in practise. Thus, either LNA input noise has to be 
significantly suppressed or noise reflections between 
the antenna and receiver have to be understood and 
parametrised with a required precision. Nevertheless, a 
better match implies easier parametrisation of the re¬ 
flections in the antenna-receiver system. Hence, it is im¬ 
portant to have these components matched as well as 
possible and thus an upgrade to a better matched an¬ 
tenna (Section 15.2|) is one of the major modifications 
needed in the presented system. 

5 Future improvements in the system 

Based on experiences with the prototype system we 
identified several weak points and used this knowledge 
to improve the system. The modifications described in 
this section have been tested in the laboratory, and a 
field deployment of the system implementing all the de¬ 
scribed changes is planned in the future. 

5.1 Modifications in the front-end 

It has been identified that reflection coefficient of the 
output (S22) of a single LNA (ZX60-33LN-S-I-) is sig¬ 
nificantly affected by the source impedance connected 
to its input (Sec. 14.31) . Therefore, in order to suppress 
effects of the source impedance changing the reflection 
coefficient of the output of the front-end (downstream 
propagation effects) more attenuators and amplifiers 
were added into a modified front-end, which introduce 
a desired “impedance isolation” between its input and 
output. Instead of having a single LNA, a front-end with 
two LNAs (ZX60-33LN-S+) providing ~44dB of gain 
separated by a 10 dB attenuator and another 10 dB at¬ 
tenuator at the output of the second LNA, which in 
fact remove all the extra added gain, was developed and 
tested in the laboratory and under the sky (but not de¬ 
ployed for longer data collection). In this configuration 
the reflection coefficient S22 measured at the output 
of the front-end is almost exactly the same whether the 
input of the new front-end is terminated, opened or con¬ 
nected to the antenna (Fig. ES|) . The initial tests of the 
new front-end in the chamber and under the sky showed 
that the Ri 1.275 MHz ripple and any other effects down 
from the output of the front-end can be calibrated out. 
Further measurements and analysis will assess whether 
suppression of impedance mismatch propagation effects 
is already sufficient in this configuration. 

Furthermore, a noise source was added to enable 
“hot/cold load” calibration schema and continuously 
determine instantaneous receiver noise temperature 
without relying on a laboratory measurement and the 
assumption that it remains constant over the mea¬ 
surement duration. Additionally, the two-position lossy 
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Figure 25. The reflection coefficient (S22) of the new front-end 
measured at the output with terminated or opened input port, 
which are the two “extremes” and the antenna is somewhere in 
between. In this case the match is very good and almost exactly 
the same in both cases. 


(~0.75 dB) solid-state RF-switch has been replaced 
with a very low loss mechanical RF-switch. A solid state 
RF-switch was used in the prototype system mainly 
for practical reasons. At the time we wanted to power 
the front-end with a single 5 V voltage (required by the 
LNA) whilst off-the-shelf mechanical RF-switches typ¬ 
ically require at least 12 V and draw much larger cur¬ 
rent. The solid-state switch, besides its relatively high 
loss, is a more complex device (possibly not entirely 
passive) than a mechanical RF-switch. Based on expe¬ 
riences with the prototype system we decided we should 
replace it by a simpler mechanical RF-switch in order to 
further simplify the signal path and remove undesired 
attenuation at the beginning of the signal path. The 
above changes should lead to significant improvements 
in the calibration procedure and be more robust to the 
effects of the source impedance propagation effects. 

5.2 Conical log spiral antenna 

The off-the-shelf biconical antenna is going to be 
replaced with a bespoke conical log spiral antenna 
(Fig. [211), which is better matched over a wide fre¬ 
quency band. Such an antenna will allow us to char¬ 
acterise all its components which should significantly 
improve the accuracy of the absolute calibration. More¬ 
over, it is a left-hand circularly polarised (LHCP) an¬ 
tenna which will help in suppressing reception of un¬ 
desired right-hand circularly polarised (RHCP) signals 
from the ORBCOMM satellites. With a narrower beam 
(FWHM Ri 70 — 80° at all frequencies) and low gain to¬ 
wards the horizon (Fig. [27l) . it will also enable us to 
use only data collected when the antenna is pointing 
towards the coldest regions of the sky. 

Initial tests showed that the antenna is well matched 
(reflection coefficient ~-20dB) in a wider frequency 
range ~ 70 — 180 MHz. Therefore, equipped with a bet¬ 
ter matched antenna and experience gained with a pro- 



Figure 26. The conical log spiral antenna (before painting) dur¬ 
ing tests at Curtin University. 


90 



Figure 27. Comparison of simulated antenna pattern of the bi¬ 
conical and conical log spiral antennas at 110, 150 and 190 MHz, 
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totype system we should be able to parameterise the 
reflected noise without any undesired attenuation be¬ 
tween the antenna and receiver. Thus, we plan to min¬ 
imise losses between the antenna and the first LNA due 
to cable, RF-switch and in particular remove the 3dB 
attenuator. 


6 Conclusions 


The precision required for the global EoR experiment is 
certainly a very challenging task from the instrumental 
and engineering points of view, for exceeding specifi¬ 
cations of typical low radio frequency applications (in¬ 
cluding radio astronomy). It requires unprecedented cal¬ 
ibration precision and thus pushes boundaries of the ra¬ 
diometer techniques at low radio frequencies. The pro¬ 
totype BIGHORNS system allowed us to learn multi¬ 
ple lessons, which triggered modifications in the initial 
signal path and will inform future deployments. The 
experiences from several deployments and laboratory 
tests once again confirmed the importance of a radio 
quiet locations for the global EoR experiment. Never¬ 
theless, even in superb radio quiet locations there are 
undesired effects such as tropospheric ducting or high 
power ORBCOMM communication satellites which sig¬ 
nificantly affect the data and require excision of a sig¬ 
nificant fraction of the data (at least 20%). In the case 
of the data presented here the excision rate was ~20%, 
but it may become even higher once more rigorous cri¬ 
teria are applied to further suppress undesired effects of 
ORBCOMM satellites transmission. 

Moreover, in the best parts of the year (when the 
Galactic Centre is in phase with the Sun) it is possi¬ 
ble to collect maximally (without taking into account 
any excision due to RFI or other effects) about 5 hours 
of data in optimal conditions (elevation of the Galac¬ 
tic Centre < —20°) and thus accuracy ^10 mK can be 
achieved in less than two nights. On the other hand in 
the worst parts of the year even 6 nights are required to 
integrate down to ~10 mK. The shorter integration time 
the better because it is much more likely to get two con¬ 
secutive nights with similar conditions (weather, iono¬ 
sphere etc.). Due to the above, the system itself must be 
extremely stable and not introduce additional variabil¬ 
ity which is already present due to environmental and 
propagation (ionosphere) factors. The stability of the 
existing system has been extensively tested and seems 
to be sufficient to integrate down the thermal noise to 
the desired precision (Figures [20] and [23|). Nevertheless, 
further sophisticated data analysis is required in order 
to separate the global EoR signal from the foregro unds 
(for example such as presented bv iLiu et al.l ( 2013!) '). 

The self-generated emission from the system was 
measured in an anechoic chamber in order to ensure 
that detection of the global EoR will not be compro¬ 
mised by the system’s own noise. Although no particu¬ 


larly significant emissions have been identified, the re¬ 
sults of the tests reassured us that these issues should 
be taken seriously in an experiment requiring such ex¬ 
treme precision. In the case of the presented system the 
^100 m separation between the antenna and the back¬ 
end electronics should be sufficient to suppress the EMI 
to a level not affecting the global EoR detection. How¬ 
ever, if possible, any amount of extra shielding (terrain 
or shielded building) should be used to further suppress 
the noise generated in the back-end electronics. 

Customised software has been developed to automati¬ 
cally reduce and calibrate the data, control their quality 
and select integrations which can be used for the global 
EoR data analysis. In order to test understanding of the 
sky signal and assess quality of the calibration, a pro¬ 
cedure integrating the FEKO antenna model with the 
sky model was developed. The results of the absolute 
calibration procedure agree with the model predictions 
to within 10%, which can be attributed to several fac¬ 
tors such as incomplete knowledge of the off-the-shelf 
biconical antenna characteristics and noise reflections 
between the antenna and the receiver, uncertainty of 
the sky model or inaccuracies in the calibration proce¬ 
dure. 

Based on the collected experiences, several weak 
points in the existing signal path, which limit calibra¬ 
tion accuracy, have been identihed. They are mainly 
due to “downstream propagation effects” related to the 
impedance of the source (antenna or 50 H terminator) 
connected to the front-end changing its output prop¬ 
erties (S22) and affecting standing noise waves down¬ 
stream in the signal path such as ripple due to the 
100 m cable. This ripple cannot be fully calibrated out 
and can only be modelled which significantly compli¬ 
cates the data analysis. Initial tests have proven that 
the problem can be dealt with by adding more com¬ 
ponents (attenuators and amplifiers) to the front-end 
which more efficiently “isolate” its output from the in¬ 
put. This modification together with several other im¬ 
provements in the front-end receiver such as change of 
solid-state RF-switch to a mechanical RF-switch and 
addition of the noise source implementing “hot/cold” 
calibration schema should significantly improve accu¬ 
racy of the calibration. The main reasons to replace 
the solid-state RF-switch with a mechanical RF-switch 
were to remove undesired attenuation at the very begin¬ 
ning of the signal path and to further simplify the signal 
path. Generally our experiences with the prototype sys¬ 
tem indicate that simplicity of the signal path is essen¬ 
tial in the global EoR experiment. Even a simple system 
has many internal and external parameters which have 
to be controlled to high precision, and introduction of 
more complex elements leads to more complex effects 
in the system. 

Finally, initial tests (such as reflection coefficient 
measurements) of the custom designed and built coni- 


PASA (2015) 

doi: 10.1017/pas. 2015.XXX 







BIGHORNS 


19 


cal log spiral antenna developed specifically for this ex¬ 
periment are very promising. It is expected to be well 
matched in the wider frequency band. Moreover, hav¬ 
ing a bespoke instrument provided more control over 
its components (such as balun, cables etc.). Therefore, 
they can be accessed, changed or characterised easier 
than in the case of the off-the-shelf product. 

The BIGHORNS system in its current form allows 
the collection of good quality data to probe the possi¬ 
bility of detecting the global EoR signal. The data col¬ 
lected in remote locations are being analysed. We expect 
that the final system equipped with a conical log spiral 
antenna and an improved receiver will sufficiently im¬ 
prove calibration accuracy to constrain the cosmological 
models or even detect the Epoch of Reionisation. 
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